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ABSTRACT 

Aims. We generate theoretical ultraviolet and extreme-ultraviolet emission line ratios for Oiv and show their strong versatility as 
electron temperature and density diagnostics for astrophysical plasmas. 

Methods. Recent fully relativistic calculations of radiative rates and electron impact excitation cross sections for O iv, supplemented 
with earlier data for A- values and proton excitation rates, are used to derive theoretical O iv line intensity ratios for a wide range of 
electron temperatures and densities. 

Results. Diagnostic line ratios involving ultraviolet or extreme-ultraviolet transitions in O iv are presented, that are applicable to a 
wide variety of astrophysical plasmas ranging from low density gaseous nebulae to the densest solar and stellar flares. Comparisons 
with observational data, where available, show good agreement between theory and experiment, providing support for the accuracy of 
the diagnostics. However, diagnostics are also presented involving lines that are blended in existing astronomical spectra, in the hope 
this might encourage further observational studies at higher spectral resolution. 
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1. Introduction 



Ultraviolet and extreme-ultraviolet emission lines arising from 
transitions in B-like O iv are detected from a wide variety of as- 
tronomical sources, ranging from the Sun (Sandlin et al. 1986) 
to other stars (Christian et al. 2004), gaseous nebulae (Feibelman 
1997) and supernova remnants (Blair et al. 1991). The diagnos- 
tic potential of these lines to provide electron temperature (T c ) 
and density (N e ) diagnostics for the emitting plasma was first 
shown by Flower & Nussbaumer (1975), who also calculated 
radiative rates plus electron and proton impact excitation cross 
sections for the ion. Since then, many authors have generated 
atomic data for O iv that have subsequently been used to derive 
theoretical diagnostic line ratios (see Tayal 2006 and references 
therein). 

Very recently, Aggarwal & Keenan (2008) have employed 
the fully relativistic grasp and Dirac rmatrx codes to calculate 
radiative rates and electron impact excitation cross sections, re- 
spectively, for all transitions among the energetically lowest 75 
fine-structure levels of O iv. These results are the most extensive 
currently available for O iv, and also should be the most reliable, 
at least for the excitation cross sections, as discussed in detail 
by Aggarwal & Keenan. In this paper we use these data, supple- 
mented with previous highly accurate calculations for radiative 
rates and proton excitation cross sections, to derive theoretical 
O iv ultraviolet and extreme-ultraviolet emission line ratios ap- 
plicable to a wide range of astrophysical plasmas. We demon- 
strate the versatility of O iv plasma diagnostics, which can pro- 
vide temperature and density estimates for a wide variety of as- 
tronomical sources ranging from low density gaseous nebulae up 
to the densest solar and stellar flares. 



2. Atomic data and theoretical line ratios 

The model ion for O iv consisted of the 75 fine-structure lev- 
els arising from the 2s 2 2p, 2s2p 2 , 2p 3 , 2s 2 3£ (£ = s, p, d), 
2s2p3£ (I = s, p, d) and 2s 2 4^ (t = s, p, d, f) configura- 
tions. Energies for all these levels were obtained from the com- 
pilation of experimental values by the National Institute of 
Standards and Technology, which may be found at their web- 
site http://physics.nist.gov/PhysRefData/ Test calculations in- 
cluding additional levels, such as those arising from the 2s 2 5^ 
and 2s2p4f configurations, where found to have a negligible ef- 
fect on the theoretical line ratios considered in the present paper. 

Einstein A-coefficients for transitions in O iv were obtained 
from the following sources: (i) Galavfs et al. (1998) for the 
forbidden 2s 2 2p 2 P[/2-2s 2 2p 2 P3/2 transition; (ii) Correge & 
Hibbert (2002) for the 2s 2 2p 2 P y -2s2p 2 4 P r intercombination 
lines; (iii) Correge & Hibbert (2004) for allowed and intercom- 
bination lines among the 2s 2 2p, 2s2p 2 , 2p 3 and 2s 2 3£ (t = s, p, 
d) levels; (iv) Aggarwal & Keenan (2008) for all remaining tran- 
sitions. For electron impact excitation rates, we have adopted 
the results of Aggarwal & Keenan, which contain several im- 
provements over previous calculations for O iv undertaken with 
the rmatrx code (Zhang et al. 1994; Tayal 2006), including the 
greatest number of levels and the largest range of partial waves. 
They are hence probably the most accurate currently available 
for this ion, as discussed in detail by Aggarwal & Keenan. 
However, there is also scope for improvement, mainly due to 
the fact that the wavefunctions adopted in the calculations of 
Aggarwal & Keenan are not as accurate as those of some other 
workers, such as Tachiev & Froese-Fischer (2000) and Correge 
& Hibbert (2002). Indeed, this is why we have adopted the A- 
values of Correge & Hibbert (2002, 2004) where possible, as 
their results have an estimated uncertainty of only +5%, com- 
pared to +20% for those of Aggarwal & Keenan. The limita- 
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tions in the wavefunctions of Aggarwal & Keenan may directly 
affect the subsequent determination of excitation rates, both for 
allowed and forbidden transitions. This is because for the weak 
allowed transitions, the A-values of Aggarwal & Keenan differ 
by up to 50% with those of Tachiev & Froese-Fischer (2000) 
and Correge & Hibbert (2004), while some of the energy lev- 
els show discrepancies of up to 8% with the experimental values 
(see Tables 1 and 2 of Aggarwal & Keenan). Nevertheless, the 
estimated accuracy of the Aggarwal & Keenan excitation rate 
data is +20% for a majority of transitions, and should be the 
most reliable currently available for O iv. 

As noted by, for example, Seaton (1964), excitation by pro- 
tons will be important for the 2s 2 2p 2 Pi/2-2s 2 2p 2 P3/2 transi- 
tion, and in the present analysis we have used the theoretical 
results of Foster et al. (1996). However, Flower & Nussbaumer 
(1975) have pointed out that proton excitation should also be in- 
cluded for the 2s2p 2 4 Py-2s2p 2 4 Py transitions in B-like ions, 
and for these rates we have adopted the calculations of Foster et 
al. (1997). Both the Foster et al. (1996) and Foster et al. (1997) 
data are estimated to be accurate to ±10%. 

Using the above atomic data in conjunction with a recently 
updated version of the statistical equilibrium code of Dufton 
(1977), relative Oiv level populations and hence emission line 
strengths were calculated for two grids of electron temperature 
(T e ) and density (N e ) values. The first grid (for T e = 10000, 
15000, 20000 and 30000 K; N e = 10 2 -10 6 ctrT 3 in steps of 0.1 
dex) is appropriate to nebular plasmas, while the second (T c = 
10 4.8_ 10 5.6 K in steps of o.i dex; N e = 10 8 -10 14 ctrT 3 in steps 
of 0.1 dex) is for solar/stellar plasma conditions. In particular, 
the adopted temperature range for the latter covers that over 
which O iv has a fractional abundance in ionization equilibrium 
of N(Oiv/N(0) > 0.006 (Bryans et al. 2008), and hence should 
be appropriate to most coronal-type plasmas. Our results are too 
extensive to reproduce here, as with 75 fine-structure levels in 
our calculations we have intensities for 2775 transitions at each 
of the 530 possible (T c , N e ) combinations considered. However, 
results involving any line pair, in either photon or energy units, 
are freely available from one of the authors (FPK) by email on 
request. Given errors in the adopted atomic data of between +5% 
and ±20% (see above), we would expect our theoretical ratios to 
be accurate to better than ±15%. 
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Fig.l. The 1(1407.3 A)/I( 140 1.1 A) line intensity ratio in Oiv, 
where I is in energy units, plotted as a function of logarithmic 
electron density (N e in ctrT 3 ) at electron temperatures of T f = 
10000, 15000, 20000 and 30000 K. 
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Fig. 2. Same as Fig. 1, but for the 1(1404.7 A)/I(1401.1 A) ratio. 



3. Results and discussion 

The Oiv intercombination multiplet around 1400 A provides 
excellent electron density diagnostics for both nebular and 
solar plasmas, as shown in Figs. 1-4, where we plot the 
emission line ratios Rj = 1(1407.3 A)/I(1401.1 A), R 2 = 
1(1404.7 A)/I(1401.1 A) and R 3 = 1(1407.3 A)/I(1404.7 A) as a 
function of T e and N e . The transitions corresponding to these 
lines are listed in Table 1, as are those for the other Oiv fea- 
tures discussed in this paper. We note that the ratios R4 = 
1(1399.7 A)/I(1401.1 A) and R 5 = 1(1397.2 A)/I(1401.1 A) have 
the same T e and N e dependence as Ri and R2, respectively (due 
to common upper levels), but with R4 = 1.02xRi and R5 = 
0.130xR 2 . 

The above theoretical ratios are in good agreement with 
observations for gaseous nebulae, such as those by Keenan et 
al. (1993). For example, for the planetary nebula NGC7662, 
Keenan et al. measured Ri = 0.30, R 2 = 0.69 and R 4 = 0.29, 
implying electron densities of log N e = 3.5, 4.0 and 3.5, respec- 
tively, from Figs. 1 and 2. These values are consistent, and also 
in agreement with the electron densities derived for NGC 7662 



Table 1. Summary of Oiv transitions. 



Wavelength (A) 


Transition 




554.51 


2s 2 2p 2 P 3/2 - 


-2s2p 2 2 P 3/2 


625.85 


2s2p 2 4 P 5/2 - 


-2p 3 4 S 3/2 


779.91 


2s2p 2 2 D 5/2 


-2p 3 2 D 5/2 


787.71 


2s 2 2p 2 P, /2 - 


-2s2p 2 2 D 3/2 


790.20 


2s 2 2p 2 P 3/2 - 


-2s2p 2 2 D 5/2 


1343.5 


2s2p 2 2 P 3/2 - 


-2p 3 2 D 5/2 


1397.2 


2s 2 2p 2 P 1/2 - 


-2s2p 2 4 P 3/2 


1399.7 


2s 2 2p 2 P 1/2 - 


-2s2p 2 4 P 1/2 


1401.1 


2s 2 2p 2 P 3/2 - 


-2s2p 2 4 P 5/2 


1404.7 


2s 2 2p 2 P 3/2 - 


-2s2p 2 4 P 3/2 


1407.3 


2s 2 2p 2 P 3/2 - 


-2s2p 2 4 P 1/2 



from line ratios in species with similar ionization potentials and 
hence spatial distributions to Oiv, such as 1(47 1 1 A)/I(4740 A) 
in Ar iv and 1(2424 A)/I(2422 A) in Ne iv, both of which indicate 
log N c = 3.5 (Keenan et al. 1997, 1998). In addition, we note 
that Keenan et al. (2002) have measured Oiv line ratios for the 
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Fig.3. The 1(1407.3 A)/l(1401.1 A) line intensity ratio in Oiv, 
where I is in energy units, plotted as a function of logarithmic 
electron density (N e in cirT 3 ) at the electron temperature of max- 
imum O iv fractional abundance in ionization equilibrium, T e = 
10 5 2 K (Bryans et al. 2008), plus +0.2 dex about this value. 
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Fig. 4. Same as Fig. 3, but for the 1(1407.3 A)/I( 1404.7 A) ratio. 



symbiotic star RRTel. The diagnostic ratios for other species 
in the RR Tel spectrum, ranging from Al n to O v, indicate log 
N e 5-8 (see Keenan et al. 2002 and references therein), over 
which density interval the theoretical values of Ri through R5 
are predicted to be effectively constant (see Figs. 1 and 2). Once 
again, there is excellent agreement between theory and observa- 
tion, with measured and predicted values of (theory in brackets) 
Ri = 0.16+0.02 (0.18), R 2 = 0.54+0.05 (0.60), R 4 = 0.17+0.02 
(0.18) and R 5 = 0.077+0.008 (0.077). 

In existing solar spectra, such as those from the HRTS and 
SOHO/SUMER instruments, the Oiv 1404.7 A line is signifi- 
cantly blended with S iv (Brage et al. 1996; Keenan et al. 2002), 
although we note that the blending is negligible under nebu- 
lar plasma conditions (Keenan et al. 2002). Additionally, the 
1397.2 A feature is very weak in solar spectra, and possibly 
blended, and hence its intensity should usually be considered an 
upper limit (Brage et al. 1996). However, the remaining useable 
density diagnostic line ratios Rj and R4 do provide consistent 
derivations of N e . For example, for Active Region B at +2" rel- 



Fig.5. Same as Fig. 3, but for the 1(1343.5 A)/I( 1407.3 A) ratio. 



ative to the limb, Feldman & Doschek (1978) measured Ri = 
0.20 and R 4 = 0. 19 from Skylab/S082A spectra, which both im- 
ply log N e = 10. 1 from Fig. 3, in good agreement with the value 
of log N e = 10.2 found from line ratios in N iv, which is formed 
at the same temperature as Oiv (Keenan et al. 1994). 

The intercombination line ratios in O iv only provide use- 
ful electron density diagnostics for values of N e up to about 
10 12 cm" 3 (see Figs. 3 and 4). However, Kastner & Bhatia (1984) 
have pointed out that the allowed lines of O iv arising from 2s2p 2 
2 P-2p 3 2 D transitions lie only ~60A from the intercombina- 
tion multiplet, and their intensity ratios are very sensitive to the 
electron density for N c > 10 11 cm" 3 . They hence should provide 
good Ne-diagnostics for very high density astronomical plasmas 
such as solar and stellar flares. In Fig. 5 we plot the ratio R$ = 
1(1343.5 A)/I(1407.3 A) as a function of both T c and N e , as the 
1343.5 A transition is the only line in the 2s2p 2 2 P-2p 3 2 D mul- 
tiplet which is unblended in existing solar spectra (Cook et al. 
1994). The measured values of R 6 = 0.89 and 0.90 for the flares 
of 1973 August 9 and 1973 September 7, from Skylab/S082B 
spectra (Cook et al. 1994), both indicate log N e = 12.5 at the 
temperature of maximum fractional abundance for O iv in ion- 
ization equilibrium, T e = 10 5 2 K (Bryans et al. 2008). These are 
consistent with the values of log N e 12.3 determined for high 
density flares using extreme-ultraviolet lines of O v, formed at a 
similar temperature to Oiv (Keenan et al. 1991). 

For temperature diagnostics, Flower & Nussbaumer (1975) 
have shown that the intensity ratio of lines within the 2s 2 2p 
2 P-2s2p 2 2 D multiplet at -790 A to those in 2s 2 2p 2 P-2s2p 2 
2 P at ~ 554 A allows T e to be estimated for the O iv emitting 
region of a plasma. This is shown in Fig. 6, where we plot 
the R 7 = 1(790.20 A)/I(554.51 A) ratio as a function of T e and 
N e . Similarly, Curdt et al. (1997) have noted that the R 8 = 
1(779.91 A)/I(787.71 A) ratio is a T c -diagnostic, and this is plot- 
ted in Fig. 7. An inspection of the two figures reveals that R 7 
is in principle a better temperature diagnostic than R 8 , as it is 
less sensitive to the adopted value of N e . However, the emission 
lines in Rg are much closer in wavelength and hence the ratio 
is more likely to be reliably measured. Indeed, Curdt et al. have 
determined R 8 = 0.033 from a SOHO/SUMER spectrum of the 
solar disk, which indicates T e ^ 10 5 1 -10 53 K from Fig. 7 (the 
exact value depending on the adopted density), in good agree- 
ment with the temperature of maximum fractional abundance in 
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Fig. 6. The 1(790.20 A)/i(554.51 A) line intensity ratio in Oiv, 
where I is in energy units, plotted as a function of logarithmic 
electron temperature (T e in K) at logarithmic electron densities 
(N e in cm" 3 ) of log N e = 8-13. 



Fig. 8. Plot of the Oiv line intensity ratio 
1(779.91 A)/I(787.71 A) against 1(625.85 A)/I(790.20 A), 
where I is in energy units, for logarithmic electron temperatures 
(T e in K) of log T e = 4.8-5.6, and logarithmic electron densities 
(N e in cm" 3 ) of log N e = 8-13. 
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Fig. 7. Same as Fig. 6, but for the 1(779.91 A)/I(787.71 A) ratio. 



ionization equilibrium for O iv, T c = 10 5 2 K (Bryans et al. 2008). 
Unfortunately, to our knowledge, the R7 ratio has not been mea- 
sured due to the low spectral resolution of existing solar ob- 
servations spanning the 554-790 A wavelength range. However, 
O'Shea et al. (1996) have determined values of the multiplet 
intensity ratio 1(790 A)/I(554 A) for several solar features from 
spectra obtained with the S-055 spectrometer on Skylab. These 
measurements lie in the range 0.47-0.56, and imply tempera- 
tures within ~ 0.2 dex of that of maximum fractional abundance 
for O iv. 

Under solar plasma conditions, O iv extreme-ultraviolet line 
ratios containing a component of the 2s2p 2 4 P-2p 3 4 S mul- 
tiplet at ~625A are strongly density sensitive (O'Shea et al. 
1996). Unfortunately, in existing solar spectra the transitions — 
which lie at wavelengths of 624.62, 625. 13 and 625.85 A — are 
all blended with the strong Mgx resonance line at 624.94 A. 
However, if the O iv features could be resolved, then in conjunc- 
tion with other relatively nearby O iv lines they would allow the 
temperature and density of the O iv emitting region of a plasma 



to be simultaneously determined via ratio - ratio diagrams, as 
illustrated in Fig. 8. 

In summary, we see that ultraviolet and extreme-ultraviolet 
emission lines of O iv provide a diverse portfolio of T e and N e 
diagnostics, applicable to a wide variety of astronomical sources 
ranging from gaseous nebulae to high electron density stellar 
flares. The present line ratio calculations, which include the most 
up-to-date atomic physics data, show good agreement with ob- 
servations where available, hence providing support for their ac- 
curacy. However, higher spectral resolution observations of sev- 
eral O iv lines, such as the components of the 2s2p 2 4 P-2p 3 4 S 
multiplet at ~ 625 A, would be very useful to allow their intensi- 
ties to be reliably measured and employed as diagnostics. 
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